Accuracy of the Hartree-Fock method for Wigner molecules at high magnetic fields 
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Few-electron systems confined in two-dimensional parabolic quantum dots at high magnetic fields 
are studied by the Hartree-Fock (HF) and exact diagonalization methods. A generalized multicenter 
Gaussian basis is proposed in the HF method. A comparison of the HF and exact results allows us 
to discuss the relevance of the symmetry of the charge density distribution for the accuracy of the 
HF method. It is shown that the energy estimates obtained with the broken-symmetry HF wave 
functions become exact in the infinite magnetic-field limit. In this limit the charge density of the 
broken-symmetry solution can be identified with the classical charge distribution. 



I. INTRODUCTION 

Properties of electron systems confined in 
quantum dots at high magnetic fields have re- 
cently become a subject of intensive theoretical 
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These studies were inspired on the one hand by the 
experimental investigation of the addition spectra of 
vertical quantum dots 2 ^ at high magnetic fields, which 
revealed a rich structure of magnetic-field induced 
ground-state transformations in the confined electron 
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and on the other hand by the search for a new symmetry 
in few-electron systems. One of the most interesting 
problems in this research is the possibility of the 
formation of Wigner moleculea 2 i 3 i 4 i 5 i 6 i 7 i 8 iV°i 25 i 26 i 27 
in which the confined electrons are distinctly spa- 
tially separated. Previous theoretical studies are 
based on the Hartree-Fock (HF) method 2 ! 3 ; 4 ! 5 ! 6 ! 8 , 
density functional theory (DFT)£, and the exact diag- 
onalization (ED) schemo 11 ! 12 . 13 ! 14 ! 15 ! 16 . 17 ! 18 ! 19 ! 20 ! 21 . 

The model confinement potential used most 
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two-dimensional (2D) cylindrically symmetric harmonic 
oscillator potential, which is a reasonable approximation 
of the confinement potential 2 ^ in vertical quantum 
dotsS. 

The external magnetic field induces ground-state 
transformations in the quantum-dot confined iV-electron 
system, which are associated with changes of the to- 
tal angular momentum and the total spin. At a cer- 
tain magnetic field the electrons become spin polarized 
and occupy the lowest-energy Fock-Darwin states with 
the z component of the angular momentum changing 
from to (1 — N)hi^. This state is called a maximum 
density droplet (MDD) 1 . In the MDD the z compo- 
nents of the total spin and total angular momentum 
take on the absolute values Nh/2 and N(N — l)h/2, 
respectively. The MDD ground state is predicted by 
the ED method 16 i 17 i 18 i 19 i 20 as well as by the HF^ and 
DFT— methods. Higher magnetic fields lead to a decay 
of the MDD, which has been observed experimentally 23 . 
This decay, considered in the framework of the ED 
methodi 6 ^ 7 -, is related to the increase of the absolute 



value of the total angular momentum above the value 
corresponding to the MDD phase. The unrestricted 
H F 2 i 3 ! 4 ! 5 and DFT£ methods predict that the ground 
state of the electron phase created after the decay of 
the MDD possesses a charge density which does not re- 
produce the symmetry of the external confinement po- 
tential. These states will be called "broken-symmetry 
states" throughout the present paper. In the broken- 
symmetry solutions obtained by the HF and DFT meth- 
ods, the electrons become localized at separate space sites 
forming a Wigner moleculo 2 ! 5 ! 6 ! 7 ! 9 ! 10 ! 25 ! 26 ! 27 . In the ED 
2 method the separation of the electrons, i.e., the forma- 
tion of the Wigner molecules, appears in the relative co- 
ordinates of the electron system and is not necessarily 
related to the broken symmetry of the charge density^. 
The HF broken-symmetry solutions are not eigenfunc- 
tions of the angular momentum operator, but are degen- 
erate with respect to rotations, i.e., can be oriented at 
an arbitrary angle. The rotational symmetry of the HF 
broken-symmetry solutions for the few-electron system 
can be restored with a post-HF treatment^. 

A multicenter basis with the one-electron wave func- 
tions was used in the theoretical studies of the 2D Wigner 
crystals 2 ^ .. Recently, a similar multicenter basis was 
use d 4 i 5 i 6 i 7 to study the quantum-dot confined electron 
systems. The paper o 4 i 5 i 6 i 7 were based on the unrestricted 
HF method 4 *^ 6 - and the variational method 7 - with the 
trial wave function in the form of a single Slater deter- 
minant with non-orthogonal one-electron wave functions. 
These one-electron wave functions 4 ^^ 7 - at high magnetic 
field yield point-like charge density distributions. There- 
fore, in the limit of infinite magnetic field, the ground- 
state charge density, obtained with the multicenter basis, 
is identical with the lowest-energy classical configuration 
of the point charges^ 1 * 3 ^ 3 .. The multicenter basis is a 
very efficient tool for the investigation of the Wigner 
molecules, since it requires only a single basis function 
per electron, while the convergence of the HF energy es- 
timates in the one-center Fock-Darwin basis^ is very slow 
in the regime of the island- like Wigner localization. Using 
the unrestricted HF method with the multicenter basis 
(MCHF), Szafran et al£ have shown that the external 
magnetic field leads to transformations of the ground- 
state symmetry of the Wigner molecules. Only in the 
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high magnetic field limit the ground-state phase (isomer 
of the Wigner molecule corresponds to the configura- 
tion of electrons, which is identical with that of a classical 
system of point charges^LSSiSi. 

In the present paper, we perform a detailed study 
of the physics behind the formation of the Wigner 
molecules. In particular, we discuss the accuracy of the 
broken symmetry solutions obtained with the multicenter 
basis in comparison with the ED results for two, three, 
and four electrons. The physical interest of this study 
relies on the investigation of the quantum systems in a 
classical localization limit. 

The paper is organized as follows: in Section II we 
briefly describe the theoretical methods, in Section III 
we present our numerical results, in Section IV - the 
discussion, and in Section V - the conclusions and the 
summary. In Appendix, we derive the wave function used 
in the present calculations. 



II. THEORY 

We consider the iV-electron system confined in the 2D 
harmonic oscillator potential with frequency loo, subject 
to the external magnetic field B oriented perpendicularly 
to the quantum dot planed. We apply the MCHF and 
ED methods. In the MCHF method we assume that 
all the electrons are spin polarized by the magnetic field 
and apply the Landau gauge, i.e., A = (—By, 0,0). We 
expand one-electron wave function ^(r) of the n-th oc- 
cupied state (/it = 1, . . . , N) 

N 
i=l 

in the basis 

^ R (r) = (a/211) 1 / 2 exp{-(a/4) (r - R) 2 + 

(i(3/2)(x-X)(y + Y)}, (2) 

where R = (X, Y), r = (x,y), and a and (3 are treated 
as nonlinear variational parameters. Function @ with 
a — (3 — m*uj c /h is the wave function of the lowest Lan- 
dau level (w c — eB/m* is the cyclotron frequency and 
m* is the electron effective mass). Moreover, function © 
with a = (2m* /h)y/cj$ + lo 2 /A and (3 = m*uj c /Ti = eB/h 
is the eigenfunction of the Fock-Darwin ground state for 
the lateral parabolic confinement p otential cen tered at 
point R with the energy equal to 1%\J w 2 + uj 2 /4 (see Ap- 
pendix). The probability density associated with wave 
function J5| is the Gaussian centered around point R. 
The centers of basis functions @ are taken from scaled 
configurations {R, ciass } of classical Wigner molecules, 
i.e., {R} = Ri = crRf ass , where i = 1,.. . ,iV 5 ' 6 . The 
scaling parameter a is the third nonlinear variational pa- 
rameter used in the present approach. In Refs. [5,6] only 
two nonlinear variational parameters were used, namely, 



the scaling parameter a and a single variational parame- 
ter a = [3. In the following, we will show that the intro- 
duction of the two independent variational parameters 
in the real and imaginary part of the exponent in Eq. J2J 
leads to a significant improvement of the variational en- 
ergy estimates at finite magnetic field. Throughout the 
present paper, the basis with the restriction a = [3 will 
be referred to as " restricted" and without it as " general- 
ized". 

In the present paper we also apply the ED method. 
The ED results, that in principle are the exact solu- 
tions of the few-electron Schrodinger equation, are used 
as reference data for the estimation of the accuracy of 
the MCHF results. In the ED calculations we use the 
symmetric gauge [A = (—By/2, Bx/2, 0)], that allows us 
to exploit the angular symmetry of the one-electron wave 
functions. The Schrodinger equation for the iV-electron 
system can be separated into the center-of-mass and 
relative-coordinate equationsi£*i2*2£. The center-of-mass 
eigenproblem possesses an analytical solution. For N — 2 
the relative-motion eigenproblem can be easily solved 
numerically with an arbitrary precision using a one- 
dimensional finite difference method. This approachi^^, 
used in the present paper for the two-electron system, is 
not applicable to the systems with a larger number of 
electrons. Thus, for N — 3 and N — 4 the ED proce- 
dure is constructed according to the configuration inter- 
action method. First, we solve the Schrodinger equation 
for a single noninteracting electron with a definite an- 
gular momentum using a finite-difference approach on 
the one-dimensional mesh with 200 points. Next, we use 
the single-electron wave functions to construct Slater de- 
terminants with the required total angular momentum 
and spin. The ^-electron Schrodinger equation is diag- 
onalized in the orthonormal basis of Slater determinants 
with proper spin-orbital symmetry and the Coulomb ma- 
trix elements are integrated numerically. The basis, i.e., 
the choice of the single-electron wave functions forming 
the Slater determinants, is optimized separately for each 
state. The calculations have been performed with a pre- 
cision better than 0.01 meV. For N = 4 this precision re- 
quires the application of the basis containing up to 2000 
Slater determinants. The ED calculations of the ground 
state for N — 3 and 4 are carried out up to 20 T with 
the maximal absolute values of the angular momentum 
equal to 18 and 34 h, respectively. 



III. RESULTS 

Let us first discuss the corrections to the MCHF en- 
ergy estimates obtained with the generalized wave func- 
tion @ . We use as an example the system of six confined 
electrons. Figure shows the difference between the en- 
ergy estimates obtained with and without the restriction 
a = (3 in wave function @. In the calculations we use 
the material parameters of GaAs, i.e. m* = 0.067 m e , 
dielectric constant e = 12.9, the effective Lande factor 
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FIG. 1: Correction to the energy estimates obtained with the 
generalized MCHF method and calculated with respect to 
the results of the restricted HF method (a = (3) for the six- 
electron system with the configuration 1-5 (solid line) and 0-6 
(dashed line) as a function of magnetic field B. The dotted 
line shows the overestimation of the energy obtained for 1-5 
configuration with the value of f3 fixed at eB/h. Inset: the 
classical 0-6 and 1-5 Wigner molecules. 



g* = —0.44, and assume the confinement energy Tiluq = 3 
meV^ 5 ". The centers of basis Q are taken from scaled 
classical configurations 1-5 and 0-6 [the phases (isomers) 
of the Wigner molecule are labelled by the numbers of 
electrons localized in the subsequent rings starting from 
the innermost one]. Phase 0-6 is the one, in which the 
Wigner molecule is created^i2ii4 after the MDD break- 
down, and possesses an intermediate character. The 1- 
5 configuration is the lowest-energy configuration of the 
classical Wigner molecule^. This is also the ground- 
state configuration of the six-electron quantum system 
at high magnetic field 5 . In the magnetic field below 4.9 
T the mutlicenter bases with both the 0-6 and 1-5 config- 
urations mimic the cylindrically symmetric MDD charge 
distribution 5 . The correction obtained with the general- 
ized basis possesses a minimum in the magnetic field in- 
ducing the MDD breakdown, for which, the application of 
the generalized basis does not improve the results. How- 
ever, the generalized basis leads to considerable correc- 
tions to the energy both in the MDD stability regime and 
after the MDD breakdown. This correction falls down to 
zero in the high-magnetic field, in which both the gen- 
eralized and restricted bases work with nearly the same 
precision. 

The dotted line in Fig. ^ shows the overestimation of 
the energy obtained for 1-5 configuration with fixed value 
of parameter [3 — eB /Ti. The value of (3 has an influence 



FIG. 2: Optimal values of the variational parameters for the 
1-5 configuration obtained with the generalized (a 7^ (3, solid 
lines) and restricted (a = /3, dashed lines) basis. Parameters 
a and (5 are expressed in units eB /h, parameter a is dimen- 
sionless. 



on the MCHF charge density due to the interference of 
the single-electron wave functions [Eqs. (|1I2|) ] centered at 
different sites. Fig. 2] shows that the variation of f3 has 
a large influence on the estimates of the energy in the 
MDD phase. On the other hand, in the Wigner molecule 
phase, the value of (3 can be safely put equal to eB/h. 

Figure [21 shows the magnetic-field dependence of the 
optimal variational parameters obtained for the 1-5 con- 
figuration in expansion The dependence of the scal- 
ing parameter a is qualitatively the same for both wave 
functions. This parameter grows rapidly when the MDD 
breaks down into the molecular phase. Generalized wave 
function © with a ^ (3 leads to the MDD decay at lower 
values of the magnetic field, which is visible in the depen- 
dence of a on the magnetic field. Parameter a tends to 
1 at high magnetic fields for which the quantum Wigner 
molecule takes the shape and size of its classical analog, 
i.e., {R} — > {R dass }. Parameters a and (3 decrease 
with the increasing magnetic field in the MDD regime 
(cf. solid lines for B <~ 5 T in Figure^- Just before 
the MDD breakdown a and (3 take on very close val- 
ues, which leads to the minimal overestimation of the 
energy obtained with restriction a = f3 at the MDD 
breakdown (cf. Fig. We have found that the in- 

crease of (3 above eB/h in the MDD regime makes the 
local maxima of the charge density 'sink' in the global 
flat maximum characteristic^ for the MDD phase. After 
the MDD breakdown parameter (3 quickly reaches eB/h,, 
i.e. takes on the value which corresponds to both the low- 
est Landau and Fock-Darwin levels [cf. the discussion of 
wave function (J2J in Sec. II]. At high magnetic fields the 
wave functions J2J centered around different sites stop to 
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overlap; in consequence, (3 stops to influence the MCHF 
charge density and takes the value eB/Ti, which ensures 
the equivalence of the centers of Landau orbitals. The 
present finding that the parameter (3 becomes equal to 
eB/Ti just after the MDD breakdown, in spite of the non- 
vanishing overlap is not evident a priori. We can give the 
following physical interpretation to this finding: the elec- 
trons in the Wigner molecule behave as if they occupied 
the independent one-particle Fock-Darwin ground-state 
orbitals, each of them centered around its own local min- 
imum of the potential energy. The optimal value of a 
rapidly grows after the MDD breakdown (cf. Fig. 2), 
which leads to the lowering of the energy obtained in the 
Wigner molecule regime (B > 5 T). The increase of a 
above eB/h enhances the electron localization and low- 
ers the electron-electron interaction energy. 
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FIG. 3: Total angular momentum (left scale) of the exact 
ground state of the two-electron system (dotted line) and 
MCHF expectation values calculated with the generalized 
(a 7^ 0, solid line) and the restricted (a = /3, dashed line) 
basis. The two curves marked by Ei nt show the expectation 
value of the electron-electron interaction energy (right scale) 
calculated with the MCHF methods. 

Figure |21 illustrates the MDD decay picture obtained 
with the ED and MCHF methods for the two-electron 
system. The dotted line shows the exact values of the 
z component of the total angular momentum as a func- 
tion of the magnetic field. The decay of the MDD is 
related with a stepwise decrease of the angular momen- 
tum from — h to — 2h, which appears for B — 5.45 T. 
The solid (dashed) line shows the expectation value of 
the total angular momentum obtained with the MCHF 
method using the generalized (restricted) basis. Figure 
3 shows that within the MDD stability regime both the 
present MCHF methods reproduce the exact value of the 
angular momentum. When the MDD breaks down, the 
MCHF expectation value of the angular momentum de- 
creases monotonically with the increasing B, in contrast 
to the exact stepwise behavior. The MDD breakdown ob- 



FIG. 4: Total angular momentum (left scale) of the exact 
ground state of the two-electron system (dotted line) and the 
MCHF expectation value calculated with the generalized basis 
(solid line) as functions of magnetic field B . The states with 
unpolarized spins are marked by "0". The dashed line shows 
the difference of the ground-state energy obtained with the 
MCHF and the exact energy calculated with the ED scheme 
(right scale). 



tained by the HF method is also related to a cusp of the 
interaction energy^ (cf. two upper curves in Fig. In 
the MDD regime the charge density distribution shrinks 
with increasing magnetic field, which results in an in- 
crease of the interaction energy. The transformation of 
the charge density from the droplet into the molecular 
phase occurs when the interaction energy exceeds some 
threshold value. The MCHF method with the generalized 
(restricted) basis leads to a MDD breakdown for B=5.55 
T (5.85 T). 

Figure 0] displays the exact ground-state angular mo- 
mentum for the two-electron system (dotted line) and the 
expectation value obtained within the HF method with 
the generalized multicenter basis (solid line) as well as 
the difference between the MCHF and the exact energy 
(dashed line) . The z component of the total spin is equal 
to fi with the exception of the low-magnetic field ground 
state and the state which appears just after the MDD 
breakdown. These two states possess zero spin and are 
labelled by " 0" in Figure 0] Notice that the expecta- 
tion value of the total angular momentum follows quite 
well the exact value. Moreover, the overestimation of 
the ground state energy within the MCHF method de- 
creases with increasing magnetic field. This decrease is 
non-monotonous, and the MCHF inaccuracy exhibit local 
minima for magnetic fields for which the exact angular 
momentum changes. 

The origin of this oscillatory behavior is explained in 
Fig. El which shows the MCHF energy estimate and 
the exact ground-state energy of the two-electron system 
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FIG. 5: Energy AE of the two-electron system calculated with 
respect to the lowest Landau level as a function of magnetic 
field B. Solid (dashed) curve show the exact (MCHF) results. 
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Same as Figure 4 but now for the three-electron 
The states with unpolarized spins are marked with 



calculated with respect to the lowest Landau level, i.e., 
AE = E-2x(hw c /2 + hg*Li B B/2). The MCHF estimate 
is a smooth function of the magnetic field, while the exact 
energy possesses cusps at the magnetic fields at which 
the ground-state angular momentum changes abruptly. 
For these magnetic fields the MCHF estimate is visibly 
closer to the exact energy value, which explains the local 
minima in Figure 4. 

A similar comparative study between the exact and the 
MCHF results has been made for the N — 3 and N = 4 
systems. For the system of three electrons the results are 
shown in Fig. 6, in which the ground states with the total 



spin equal to h/2 are labelled by "1/2" . The other states 
are fully spin polarized. Similarly as in the case of the 
two-electron system the ground state which appears after 
the MDD breakdown is not spin polarized. The overes- 
timation of the total energy obtained with the MCHF 
method exhibits a similar qualitative dependence on the 
magnetic field as for two electrons. In contrast to the 
two-electron case, the MCHF method with the general- 
ized basis predicts a breakdown of the MDD for a slightly 
smaller magnetic field value (B = 4.6 T) than the exact 
result (B = 4.8 T). The MCHF with the restricted basis 
yields the magnetic field inducing the MDD breakdown 
5 = 5 T. 
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FIG. 7: Same as Figure 4 but now for the four-electron sys- 
tem. The states with unpolarized spins are marked with the 
value of the z-component of the total spin in Ti units. 

A similar result for the four-electron quantum dot is 
shown in Figure[7| The non-fully-polarized ground states 
are marked by the quantum numbers of the total-spin z 
component "0" and "1". Contrary to the two and three 
electron systems, the four-electron MDD decays into a 
spin-polarized state, but a low-spin state still appears at 
the higher magnetic field. The transitional appearance of 
low-spin ground states at magnetic fields after the MDD 
breakdown have been reported first in Ref^ 1 for N = 2 
and N = 4. The exact magnetic field inducing the MDD 
decay is equal to 4.75 T, while the MCHF with the gener- 
alized basis predicts a value of 4.38 T and the MCHF with 
the restricted basis gives 4.8 T. In this case the MCHF 
with the restricted basis gives accidentally a better es- 
timate for the MDD breakdown field. For five electrons 
the difference between both the MCHF estimates of the 
magnetic field inducing the MDD decay is 0.4 T. For six 
and more electrons the differences are not larger than 0.2 
T. 

Figures 4, 6, and 7 show that the MCHF inaccuracy 
decreases with the increasing magnetic field. In order 
to find the high-field asymptotic behavior of the MCHF 
energy estimate we have plotted in Figure [S] the MCHF 
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FIG. 8: The MCHF inaccuracy for the ground state for the 
two-electron system as a function of 1/B for N=2 (solid line), 
3 (dashed line) and 4 (dash-dotted line) The dotted line shows 
the high-magnetic field asymptote for N = 2, parametrized 
as 4.026/B [meV/T]. 



error as a function of 1/B for N = 2, 3, and 4. The plot 
for two electrons covers the magnetic fields up to 160 T, 
while the plots for three and four electrons are drawn up 
to 20 T only. The results for N = 2 show that at high 
magnetic fields the MCHF inaccuracy is proportional to 
1/B. At high magnetic fields the asymptotic behavior 
of the MCHF inaccuracy for N = 2 can be very well 
approximated by the function f(B) = 4.026/5 [meV/T], 
and consequently the MCHF approach becomes exact for 
B — > oo. The plots for N = 3 and 4 in the studied 
(narrower) range of the magnetic field exhibit a similar 
tendency as that for N = 2, however, they do not become 
linear functions of 1/B for B < 20 T. The comparison of 
the MCHF inaccuracies for N — 2, 3, and 4 indicates that 
the overestimation of the MCHF ground-state energy in 
the magnetic field range above the MDD instability does 
not substantially increase with TV. For B — 20 T the 
overestimation of the ground-state energy for N = 2,3, 
and 4 is equal to 0.21, 0.24, and 0.18 meV, respectively. 



IV. DISCUSSION 

Figs. 01 El and [7\ show that at high magnetic field the 
exact ground-state angular momentum take the so-called 
magic values^Sii&Si and change by Nh. Only for these 
magic values of the angular momenta the classical sym- 
metry can be reproduced in the inner coordinates of the 
quantum systems^. On the other hand the classical sym- 
metry is ensured in the MCHF by the present choice of 
centers of orbitals @ and the linear change of the expec- 
tation value of the total angular momentum is related 
with the growing localization of wave functions |J2J . 



The results presented in Figs. 211^113 and [5] show that 
the broken-symmetry solutions obtained with the MCHF 
method provide exact energy results in the high magnetic 
field limit. This fact might be rather surprising, since one 
could expect that the exact solutions of the few-electron 
Schrodinger equation should also be the eigenfunctions of 
the angular momentum operator. The exact eigenfunc- 
tions yield the charge density distribution, which repro- 
duces the symmetry of the confinement potential. This 
apparent contradiction can be solved if we consider the 
Schrodinger equation for the electron system confined in 
the parabolic potential. This equation can be separated 
in the center-of-mass and relative coordinates. In the 
framework of the ED approach, the separation of the 
electrons, i.e., the Wigner localization, appears in the 
relative (inner) coordinates of the electron system, while 
the charge density in the laboratory frame is affected 
by the center-of-mass motion. The center-of-mass eigen- 
problem of the few-electron system has the same form 
as the Fock-Darwin equation for the single electron. In 
the high magnetic field limit, this equation possesses a 
degenerate ground state (the lowest Landau state), for 
which the angular momentum is arbitrary. In the case 
of this degeneracy a superposition of ground states with 
different angular momenta is still the ground state of the 
Fock-Darwin equation, even though it leads to solutions 
with the broken symmetry of the charge density distri- 
bution in the laboratory frame. In the HF method, the 
formation of the Wigner molecule, i.e., the separation of 
the electrons, is only possible in the broken-symmetry so- 
lutions. The results of the present paper show that exact 
energy is obtained for these broken-symmetry solutions 
in the infinite magnetic field limit. 

In the recent paper of Bcdnarek et al& a study of the 
accuracy of the HF method has been presented for the 
quasi-one-dimensional (ID) structure. In the quasi-lD 
structures the unrestricted HF method becomes exac1i2£ 
in the large quantum dots, in which the Wigner molecules 
are formed^. In the ID structures, the HF method con- 
serves the two-particle parity symmetry of the exact so- 
lution; so, there is no problem with the broken symmetry 
of the HF solutions like in the 2D circularly symmetric 
quantum dots. 

The present results show that the overestimation of the 
exact energy obtained in the broken-symmetry MCHF 
solutions is relatively small at the magnetic fields, for 
which the ground state is degenerate. The ground state 
of the few-electron system is twofold degenerate at these 
magnetic fields, which induce a stepwise change of the 
ground-state angular momentum (cf. Figs. and[7J|. 
In this case, the exact ground state can be a superpo- 
sition of two states with different angular momenta and 
therefore the charge density can have the symmetry dif- 
ferent from that of the external potential. For these fields 
the broken symmetry of the HF solutions leads to the de- 
crease of the energy separation between the MCHF and 
ED results (cf. Fig. 6), which in turn causes the charac- 
teristic oscillations of the MCHF error as a function of 
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the magnetic field as shown in Figs. 14171 

The application of the generalized variational basis re- 
sults in a modification of the phase diagram^, for the 
Wigner molecules. This modification is due to the dif- 
ferent precision of the restricted wave function for dif- 
ferent phases (cf. Fig. 1). The critical magnetic fields 
for the MDD breakdown, obtained with the generalized 
basis, are shifted toward lower values (cf. Fig. 3) and 
the range of the stability of different phases is modi- 
fied. The improved results conserve the characteristic 
features of the original phase diagram^ i.e., the interme- 
diate phases correspond to the configurations, for which a 
larger number of electrons is gathered on the outer ring 
of the molecule in comparison with the classical, high- 
magnetic-field, ground-state configuration. 

The present calculations performed for small num- 
ber of electrons indicate that in the high magnetic field 
the inner charge distribution of electrons can be derived 
from classical calculations for particles interacting via a 
Coulomb (1/r) interaction. On the other hand in the 
Laughlin wave function^, which seems to become an ex- 
act description of the many-body state at high magnetic 
field, the distribution of electrons in the inner coordinates 
corresponds to classical configuration of particles inter- 
acting with a logarithmic potential. It is therefore not 
excluded that for larger number of electrons the localiza- 
tion may be different than for electrostatically interacting 
classical particles. 



V. CONCLUSIONS AND SUMMARY 

We have investigated the quantum-dot confined 
TV-electron system at high magnetic fields using the 
HF method with the generalized multicenter basis and 
the exact diagonalization method. We have indicated 
that the magnetic-field dependence of the variational 
parameters of the generalized basis can be used as one 
of the signatures of the liquid-solid phase transition, i.e., 
the breakdown of the MDD into the molecular phase. 
The occurrence of the cusp of the interaction energy as 
a function of the magnetic field and the decrease of the 
MCHF angular-momentum expectation value below that 
corresponding to the MDD are the other signatures of 
the MDD decay. We have discussed the accuracy of the 
energy estimates obtained with the broken-symmetry 
HF solutions, which - in contrast to the exact solutions 
of the Schrodigner equation - are not eigenfunctions of 
the total angular-momentum operator. It turns out that 
the angular-momentum eigenvalues in the MDD phase 
are reproduced with the high accuracy by the MCHF 
expectation values, which at higher magnetic fields 
linearly decrease with increasing B in contrast to the 
exact stepwise decrease. The results of the present paper 
show that the MCHF inaccuracy decreases with the 
increasing magnetic field and that the MCHF method 
basis yields the exact ground-state energy in the infinite 
magnetic- field limit. We have found the characteristic 



oscillations of the HF inaccuracy, which exhibits local 
minima at those magnetic fields, for which the exact 
ground state is degenerate with respect to the angular 
momentum. The relation of these oscillations with the 
existence of the exact ground states with the broken 
symmetry has been pointed out. The results of the 
present paper show that the envelope of the MCHF 
inaccuracy oscillations at high magnetic field decreases 
linearly to as a function of 1 /B and that in the strictly 
infinite magnetic field limit the exact energy is obtained 
for broken-symmetry HF solution with the classical 
point-charge distribution. 
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Appendix 

In order to derive wave function [Eq. (J2J] let us first 
consider the single electron in a homogeneous magnetic 
field. In the Landau gauge, i.e., for A(r) = (—By, 0,0), 
the Hamiltonian has the form 

where lo c = eB/m*. The lowest Landau energy level 
E = hu c /2 is infinitely- fold degenerate. The corre- 
sponding degenerate eigenstates have the form 

Xg (x, y) = d exp[iqx - {fl/2){y - q/(3) 2 } , (4) 

where [3 — m*uj c /h — eB/h, C\ is the normalization 
constant, and q 6 (— oo,+cx)). Due to this degeneracy, 
an arbitrary linear combination of wave functions J3J 1S 
an eigenfunction of Hamiltonian J3J to eigenvalue Eq. 
The most general form of this linear combination can be 
written as 

/+oo 
f{q)xq(x,y)dq ■ (5) 
-00 

Taking 

f(q) = cxp [q(y - ix Q ) - q 2 /20\ (6) 

and performing the integration in Eq. JSJ, we obtain 

ip (x, y) = C exp{-(^/4) [(x - x ) 2 + (y - y ) 2 ] + 

(i(3/2)(x-x )(y + y )} , (7) 

where Co is the normalization constant. Wave function 
@ corresponds to the electron probability density, which 
is localized around center r = (xo 7 yo). For arbitrary r 
wave function J7J is the eigenfunction of the Hamiltonian 
© associated with eigenvalue Eq. Due to the arbitrary 
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choice of ro, wave functions J7J), localized at different 
centers, correspond to the same lowest Landau level. 

If, in addition to the magnetic field, the electron is 
subject to the external parabolic potential centered at 
site R = [X,Y), i.e., 

Vconf(x, y) = ^[(x Xf + (y- Y) 2 ] , (8) 
we deal with the Fock-Darwin eigenproblem. Then, the 



ground-state wave function is centered around R and is 
written down in the normalized form as 

Ax.Y) (x, y) = cxp{-(a/4) [(a; - Xf + (y - Y) 2 ] + 

(iP/2)(x-X)(y + Y)}/(a/2^/ 2 , (9) 



where a — (2m* /U)^J(Jq + u>% /4. Therefore, we obtain 
the wave function of form 
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